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ABSTRACT 

We present seven epochs of spectropolarimetry of the Type lib supernova (SN) 2011dh 
in M51, spanning 86 days of its evolution. The first epoch was obtained 9 days after 
the explosion, when the photosphere was still in the depleted hydrogen layer of the 
stripped-envelope progenitor. Continuum polarization is securely detected at the level 
of P « 0.5% through day 14 and appears to diminish by day 30, which is different from 
the prevailing trends suggested by studies of other core-collapse SNe. Time-variable 
modulations in P and position angle are detected across P-Cygni line features. Ha 
and He i polarization peak after 30 days and exhibit position angles roughly aligned 
with the earlier continuum, while Oi and Can appear to be geometrically distinct. 
We discuss several possibilities to explain the evolution of the continuum and line 
polarization, including the potential effects of a tidally deformed progenitor star, as- 
pherical radioactive heating by fast-rising plumes of ®®Ni from the core, oblique shock 
breakout, or scattering by circumstellar material. While these possibilities are plausi¬ 
ble and guided by theoretical expectations, they are not unique solutions to the data. 
The construction of more detailed hydrodynamic and radiative-transfer models that 
incorporate complex aspherical geometries will be required to further elucidate the 
nature of the polarized radiation from SN 2011dh and other Type Hb supernovae. 
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1 INTRODUCTION 

Spectropolarimetric observations of supernovae (SNe) probe 
the explosion geometry and the relative distribution of 
chemical elements within the progenitor’s outer envelope 
and inner ejecta (McCall 1984; see also the review by Wang 
& Wheeler 2008). This is valuable information that can pro¬ 
vide clues to the nature of the progenitor systems, the explo¬ 
sion mechanism, and the kicks imparted to newborn neutron 
stars. 

The most common intrinsic source of linearly polar- 
* E-mail: mauerhan@astro.berkeley.edu 


ized emission in SNe is Thomson scattering of photons by 
free electrons in the dense ionized outflow or circumstel¬ 
lar medium. Optical scattering of photons by circumstel¬ 
lar dust particles is another potentially important mecha¬ 
nism for efficiently producing linear polarization (Wang & 
Wheeler 1996). However, in the idealised case of spherically 
symmetric geometry, every electric-held vector from the cir¬ 
cular scattering surface, as observed on the sky, will have an 
orthogonally oriented vector of equal magnitude in an ad¬ 
jacent quadrant of the circle, and will thus cancel out. The 
detection of net polarization therefore requires some degree 
of asphericity on the plane of the sky to break the symmetry. 
Net continuum polarization indicates the presence of global 
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asphericity for the electron-scattering photosphere and/or 
circumstellar medium, or a patchy distribution of contin¬ 
uum opacity. Polarization changes across line features, on 
the other hand, can occur for both globally spherical and 
aspherical configurations, when there is partial obscuration 
of the underlying photosphere by an uneven or clumpy dis¬ 
tribution of line opacity (Kasen et al. 2003; Dessart & Hillier 
2011). These general features are potentially important for 
both thermonuclear and core-collapse (CC) explosions. How¬ 
ever, in order for spectropolarimetry to provide information 
on the intrinsic asphericity of SNe, one must have knowledge 
of the wavelength-dependent linear polarization produced by 
the differential absorption of light by interstellar dust along 
the line of sight (Serkowski 1968), from both the Milky Way 
and the host galaxy. 

A large fraction of SNe that have been observed with 
polarimetry exhibit evidence for some degree of aspheric¬ 
ity, but the spectropolarimetric properties are wide ranging. 
As studies have progressed in the last decade, systematic 
trends have emerged for the various SN subtypes. Ther¬ 
monuclear SNe la generally exhibit a relatively low degree 
of intrinsic net polarization in the continuum (P < 0.2%; 
Wang & Wheeler 2008; Wang et al. 2003a; Kasen et al. 
2003) and more substantial polarization across P-Cygni ab¬ 
sorption features (occasionally ~ 1% or more; e.g., Leonard 
et al. 2005; Patat et al. 2009), indicating only modest as¬ 
phericity for the overall geometry of the outflow, yet a some¬ 
times highly nonuniform distribution of line optical depth. 
CC SNe have exhibited markedly different evolution. Some 
of the most well-studied SNell-P have exhibited an evolv¬ 
ing degree of continuum polarization, beginning at low or 
modest levels during the plateau phase and then sharply in¬ 
creasing as the photosphere recedes below the H envelope 
and into the more highly aspherical He layer and metal-rich 
ejecta (Leonard et al. 2006; Chornock et al. 2010; but also 
see Leonard et al. 2012, and preliminary results of Leonard 
et al. 2013). SNelln, which are explosions that interact with 
circumstellar material (CSM) to produce their namesake rel¬ 
atively narrow emission lines, commonly show a moderately 
high degree of continuum polarization at early times (1-2%; 
Leonard et al. 2000; Hoffman et al. 2008; Smith et al. 2008; 
Patat et al. 2011; Mauerhan et al. 2014), and in some cases 
exhibit a complex temporal evolution in polarization magni¬ 
tude and position angle on the sky (Mauerhan et al. 2014). 
These characteristics have been interpreted as evidence for a 
superposition of multiple, geometrically distinct components 
of polarization from the SN outflow and a highly aspherical 
distribution of shocked CSM, consistent with toroidal geom¬ 
etry in several cases (Hoffman et al. 2008; Mauerhan et al. 
2014). 

Stripped-envelope SNe are explosions that stem from 
stars that have depleted hydrogen envelopes (Type Ilb) or 
that are entirely devoid of hydrogen (Type Ib) and helium 
(Type Ic) lines; see Filippenko (1997) for a review of SN 
classification. These subclasses of SNe consistently exhibit 
substantial polarization in the continuum and enhanced po¬ 
larization across lines (Leonard et al. 2002; Wang et al. 
2003b; Wang & Wheeler 2008). SNellb, in particular, are 
an intermediate class between SNell-P and Ib, and consti¬ 
tute ~ 11% of CC SNe in large Milky Way-like galaxies 
(Smith et al. 2011; Li et al. 2011). In the majority of cases, 
the stellar progenitors at the time of explosion are thought 


to have been yellow supergiants or Wolf-Rayet stars of rela¬ 
tively modest initial mass (10-20 Mq) that had the majority 
of their hydrogen envelopes stripped down to Mh ~ 0.001- 
1 Mq via gravitational interaction with a binary companion 
(e.g., Podsiadlowski et al. 1993; Woosley et al. 1994; Utrobin 
1994; Hoflich 1995; Maund et al. 2004; Claeys et al. 2011; 
Smith et al. 2011; Van Dyk et al. 2013; Gal-Yam et al. 2014; 
Jerkstrand et al. 2015), in some cases potentially forming 
stripped-envelope Wolf-Rayet stars before exploding. These 
hypotheses are supported by the direct detection of the pro¬ 
genitor stars and their candidate companions in respective 
pre- and post-explosion images from the Hubble Space Tele¬ 
scope (HST) (Maund et al. 2004; Van Dyk et al. 2011, 2013; 
Fox et al. 2014). 

Chevalier & Soderberg (2010) suggested that SNellb 
be subcategorized into two types of explosions that stem 
from either (i) progenitors with extended atmospheres and 
slow dense winds (ellb; R « 100 Rq), leading to a prominent 
early peak in the light curve from the extended shock-heated 
envelope, or (ii) more compact progenitors with fast tenu¬ 
ous winds (cllb; i? « 1 Rq). The latter subclass appears to 
be distinguishable by weak optical emission from the shock- 
heated envelope at early times, nonexistent or very weak H 
emission in the late nebular phase, rapid radio evolution, and 
nonthermal X-ray emission. The cllb class is thus more con¬ 
sistent with Wolf-Rayet-like progenitor stars than cool su¬ 
pergiants having extended atmospheres. The classification of 
SNellb within these subcategories is not always clear, how¬ 
ever, as there could be a continuum of intermediate objects 
having progenitor radii between the ellb and cllb extremes 
(Horesh et al. 2013). Different amounts of ultraviolet (UV) 
excess might also be an important diagnostic for determin¬ 
ing the nature of SNIIb progenitors and the properties of 
their CSM (Ben-Ami et al. 2015). 

At early phases, SNellb display Balmer-dominated 
spectra and commonly exhibit a low to moderate-level con¬ 
tinuum polarization (see the review by Wang & Wheeler 
et al. 2008, and references therein), indicative of modestly 
aspherical global outflow. By the time of maximum light, 
just weeks after explosion, the photospheres of SNellb have 
already passed through the depleted hydrogen layer, giving 
rise to He-dominated spectra; this Hb—J-Ib transition is often 
accompanied by a substantial increase in line polarization. 
Well-studied examples for which quality spectropolarimetry 
has been obtained include the SN Hb prototype SN 1993J 
(Trammel et al. 1993; Hoflich 1995; Hoflich et al. 1996; Tran 
et al. 1997) and the remarkably similar SN 1996cb (Wang 
et al. 2001), both of which exhibited substantial continuum 
polarizations of ~ 0.6% before peak brightness and H-rich 
phases that lasted for 2-3 weeks after the explosion. The 
spectropolarimetric similarity between these SNe is note¬ 
worthy, given that SN 1993J has been classified as an ellb 
and SN 1996cb as a cllb (Chevalier & Soderberg 2010). The 
additional case of the Type Hb SN 2001ig began with a low 
(~ 0.2%) continuum polarization at early times, and later 
rose to ~ 1% after the photosphere receded into the He layer 
(Maund et al. 2007a), similar in behaviour to SNell-P. The 
more recent cllb SN 2008ax was a rather unusual object that 
exhibited the Balmer-dominated spectrum of a young SN Hb 
only during the first several days after exploding, indicative 
of an H envelope mass that was substantially lower than 
that of most other SNellb (e.g., see Marion et al. 2014). 
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The object exhibited substantial continuum polarization of 
~ 0.6% at early times, indicative of an aspherical outer en¬ 
velope, and developed very strong Ha line polarization of 
~ 3.4%, stronger than ever observed in any SNIIb and with 
large radial velocities reaching 20,000 kms“^ (Chornock et 
al. 2011). 

Overall, spectropolarimetric properties of SNellb are 
broadly consistent with a moderately aspherical photo¬ 
sphere, punctuated by irregular clumps or plumes of en¬ 
hanced line optical depth in the ejecta. Optical scattering 
of light from CSM dust particles is also a potentially impor¬ 
tant factor affecting the degree of continuum polarization, 
and was possibly influential in the case of the prototypical 
SN1993J (Tran et al. 1997). 

SN2011dh was a bright Type lib explosion in the 
nearby spiral galaxy NGC5194 (Messier 51a, hereafter re¬ 
ferred to as M51), discovered independently by several am¬ 
ateur astronomers (summarized by Griga et al. 2011 and 
Arcavi et al. 2011). Remarkably, a yellow supergiant (YSG) 
star was found to be spatially coincident with the SN in HST 
archival images of the galaxy (Li & Filippenko 2011; Maund 
et al. 2011; Van Dyk et al. 2011). It was unclear whether this 
star was the actual progenitor of the explosion, a neighbor, 
or even a surviving companion star to a more compact pro¬ 
genitor as initially hypothesized on the basis of its early-time 
luminosity (Arcavi et al. 2011) and radio emission (Soder- 
berg et al. 2012). Eventually, the identification of the YSG 
as the true progenitor was solidified after the SN had faded, 
revealing that the YSG star had disappeared (Van Dyk et 
al. 2013). The same star was found to be variable in ground- 
based images up to 3 yr before the explosion (Szczygiel et al. 
2012). Recently, a candidate hot companion star has been 
detected in late-time UV data from HST (Folatelli et al. 
2014); SN2011dh is thus the second SNIIb to exhibit direct 
evidence for a binary companion, after SN 1993J (Maund et 
al. 2004; Fox et al. 2014). 

Estimates of the YSG progenitor’s initial main-sequence 
mass and radius lie in the range 10-19 M© and 200-300 Rq, 
respectively (Sahu, Anupama, & Ghakradhari 2013; Van 
Dyk et al. 2011; Bersten et al. 2012; Benvenuto et al. 2013; 
Ergon et al. 2014); SN2011dh is thus a member of the ellb 
subclass. The ellb designation is also consistent with the 
photometric detection of the shock-heated envelope phase in 
the early-time light curve (Arcavi et al. 2011), although the 
fact that this phase was relatively short lived compared to 
that of SN 1993J suggests that the progenitor of SN 2011dh 
was substantially less extended. Spectroscopic evolution of 
SN 2011dh showed an early H-rich phase lasting for roughly 
two weeks, after which the photosphere entered the He layer 
(Arcavi et al. 2011; Marion et al. 2014; Ergon et al. 2014). 
The late-time decline in the light curve showed the explo¬ 
sion to be consistent with the production of 0.09 Mq of ®®Ni 
and 0.2 Mq of oxygen, and also provided evidence for possi¬ 
ble late-time dust formation, GSM interaction, and substan¬ 
tial asphericities in the explosion ejecta (Sahu et al. 2013; 
Shivvers et al. 2013). Moreover, recent high-resolution VLBI 
imaging at 8.4 GHz might have detected “hot spots” poten¬ 
tially indicative of modestly aspherical distribution of the 
outer ejecta, although this feature might just be an artifact 
of the data-reduction process (de Witt et al. 2015). Spec¬ 
tropolarimetry of SN 2011dh is thus well motivated. 

Here we present high-quality multi-epoch spectropo¬ 


larimetry of SN 2011dh, producing the most complete 
dataset of its kind for a SNIIb. In §2 the observational fa¬ 
cilities are described. We present the data and analyze the 
results in §3. Section 4 compares SN 201 Idh to other SNellb 
that have been the subject of spectropolarimetric observa¬ 
tions and interprets the data in the context of existing ex¬ 
plosion and progenitor models. 

2 OBSERVATIONS 

2.1 The Kast Spectrograph at Lick Observatory 

The Shane 3 m reflector at Lick Observatory equipped with 
the Kast spectrograph (Miller & Stone 1993) was used to 
observe SN2011dh on 2011 June 6 and 29, and on Au¬ 
gust 4 and 25 (UT dates are used throughout this paper). 
Kast is a dual-beam spectropolarimeter that utilizes a ro¬ 
tatable semiachromatic half-waveplate to modulate the in¬ 
cident polarization and a Wollaston prism in the collimated 
beam to separate the two orthogonally polarized spectra 
onto the CGD detector. Only the red channel of Kast was 
used for spectropolarimetry; a GG455 order-blocking filter 
suppressed all second-order light at wavelengths shortward 
of 9000 A. Observations were made with the 300 line mm“^ 
grating and the 3”-wide slit, yielding a spectral resolution of 
~ 4.6 Apixel”^ and a full width at half-maximum intensity 
(FWHM) spectral resolution of ~ 16 A (as determined by 
the FWHM of lamp lines). The useful wavelength coverage 
with this setup is 4600-9000 A. The orientation of the slit 
on the sky was always set to a position angle of 180° (i.e., 
aligned north-south). Exposures of 900 s were obtained at 
each of four waveplate positions (0?0, 45?0, 22?5, and 67?5). 
In most cases, three waveplate sequences were performed 
and the results median combined. Flatfield and arc lamp 
(for wavelength calibration) spectra were obtained immedi¬ 
ately after each sequence, without moving the telescope. The 
arc lamps are internal to the Kast instrument, while the flat- 
held spectra were produced by the reflection of incandescent 
lamp light off of the inner surface of the dome, which was 
rotated into the field of view of the telescope before or af¬ 
ter the science exposures. The flatfield spectral images were 
normalized by a low-order spline fit to the continuum before 
dividing them into the science images. 

Eor polarimetric calibrations, standard stars were se¬ 
lected from the sample of Schmidt et al. (1992a,b). The 
unpolarized standard stars BD-|-26°2606 and BD-|-32°3739 
were observed to verify the low instrumental polarization 
of the Kast spectrograph. We constrained the average frac¬ 
tional Stokes Q and U values to < 0.05%. By observing the 
above unpolarized standard stars through a 100% polariz¬ 
ing filter, we determined that the polarimetric response is 
so close to 100% that no correction was necessary, and we 
obtained the instrumental polarization position-angle curve, 
which we used to correct the data. We observed the high- 
polarization stars HD 154445, HD 127769, HD 161056, and 
HD 204827 to obtain the zeropoint of the polarization posi¬ 
tion angle on the sky (0) and to determine the accuracy of 
polarimetric measurements, which were generally consistent 
with previously published values to within AP < 0.05% and 
AO < 1°. 
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Table 1. Integrated V-hand and continuum polarization of SN 2011dh 


UT Date 

MJD 

day'^ 

Tel./Instr. 

Pv {%)'• 

Sv (deg)*’ 

Pcont (%r 

^cont (deg)^ 

Jun. 09 

5721.8 

9 

Lick/Kast 

0.406 (0.010) 

34.9 (0.8) 

0.47 (0.02) 

23.5 (0.8) 

Jun. 14° 

5727 

14 

Bok/SPOL 

0.237 (0.013) 

14.8 ^.0) 

0.45 (0.02) 

2.3 (1.3) 

Jun. 28 

5742.8 

30 

Lick/Kast 

0.358 (0.013) 

8.7 (1.0) 

0.18 (0.04) 

12.0 (3.4) 

Jul. 06 

5748.7 

36 

P200/DBSpec 

0.345 (0.021) 

18.7 (1.7) 

0.21 (0.03) 

29.1 (2.7) 

Jul. 28“* 

5771 

58 

Bok/SPOL 

0.425 (0.042) 

18.3 (4.2) 

- 

- 

Aug. 04 

5777.7 

65 

Lick/Kast 

0.400 (0.022) 

13.2 (1.7) 

- 

- 

Aug. 25 

5798.7 

86 

Lick/Kast 

0.259 (0.035) 

15.2 (3.9) 

- 

- 


^Day is with respect to the adopted explosion date (JD 2455713.08), and ronnded to the nearest integer value. 

°14-band values averaged over the wavelength range 5050—5950 A. 

Average of data from June 14 and June 15. 

Average of data from July 28 and July 30. 

^Average of data from Aug 4 and Aug 5. 

^Continuum sample region is 7000-7450 A; some weak line contamination might be present for the June 28 and July 6 data. 


2.2 SPOL at Arizona Observatories 

Spectropolarimetric data were obtained on 2011 June 
14, 15, 28, and 30 with the Bok 90 inch telescope on 
Kitt Peak. Observations made use of the CCD Imag- 
ing/Spectropolarimeter (SPOL; Schmidt et al. 1992a). Like 
Kast, SPOL contains a rotatable semiachromatic half- 
waveplate that is used to modulate incident polarization and 
a Wollaston prism in the collimated beam to separate the 
two orthogonally polarized spectra onto the detector. SPOL 
was configured with a 600 line mm“^ grating in first order. 
The slit selection was based on the seeing conditions and 
ranged from 1.1” to 5.1” x 51”. This setup yielded a spectral 
resolution of 20-30 A (FWHM) with useful wavelength cov¬ 
erage in the range 4000-7600 A. A standard Hoya L38 block¬ 
ing filter was used to ensure that the first-order spectrum 
was not contaminated by second-order light for A > 7600 A. 
The efficiency of the waveplate as a function of wavelength 
was measured and corrected at all epochs by inserting a fully 
polarizing Nicol prism into the beam above the slit. The ori¬ 
entation of the slit on the sky was always set to a position 
angle of 0° (i.e., aligned north-south). A series of four sepa¬ 
rate exposures that sample 16 orientations of the waveplate 
yields two independent, background-subtracted measures of 
each of the normalized linear Stokes parameters, q and u. 
For some observing runs, several polarization sequences of 
SN 2011dh from different nights were combined if no po- 
larimetric change was discernible, with the weighting of the 
individual measurements based on photon statistics. 

The instrumental polarization of Bok/SPOL is very 
small (< 0.1%), as measurements of unpolarized standard 
stars over the past two decades have consistently shown. 
Still, to confirm, we observed one or both of the unpolarized 
standard stars BD-|-28°4211 and G191B2B (Schmidt et al. 
1992b) during each epoch, confirming that the instrumental 
polarization is < 0.1% and that there is no low-level residual 
structure observed in the q and u Stokes spectra. The linear 
polarization position angle on the sky {9) was determined by 
observing one or more of the interstellar polarization stan¬ 
dards Hiltner 960, VI Cyg #12, or BD+59°389 (Schmidt et 
al. 1992b) during all epochs. The adopted correction from 
the instrumental to the standard equatorial frame for 9 for 
all epochs was determined from the average position angle 
offset of the polarized standards. Differences between the 
measured and expected polarization position angles were 
< 0?3 for all of the standard stars. 


2.3 The Double Spectrograph at Palomar 
Observatory 

With the 5 m Hale telescope at Palomar Observatory, we ob¬ 
tained a single epoch of spectropolarimetry of SN 201 Idh on 
2011 July 6, using the Double Spectrograph (DBSP; Oke & 
Gunn 1982) with polarimeter. The total exposure time was 
5760 s, spanning four complete sets of spectropolarimetric 
data (i.e., each set consisted of four 360 s observations at 
each of four waveplate position angles). 

Owing to spurious polarization features being intro¬ 
duced by the dichroic (see, e.g., Ogle et al. 1999), we ob¬ 
served with only the red camera of the spectrograph (i.e., 
no dichroic) using the 158/7500 grating and the GG445 
order-blocking filter (to mitigate second-order light contam¬ 
ination). This resulted in a final, useful wavelength range of 
5100-9100 A, with a spectral resolution of ~ 15 A (FWHM). 
We note that the absence of a dichroic resulted in poor spa¬ 
tial and spectral focus in one of the beams (the “top” beam 
of the two beams delivered by the spectropolarimeter on 
the CCD chip) blueward of 7500 A, with the resolution ap¬ 
proaching ~ 30 A near the blue end of the spectral range. In 
order to preserve the best spectral resolution, the total-flux 
spectrum was produced using only the bottom beam in the 
wavelength range 5100-7500 A, with the two beams being 
combined redward of 7500 A. For the spectropolarimetry, 
we retain the entire spectral range, and so the resolution of 
any spectropolarimetric features varies somewhat from the 
blue end to the red end of the spectrum. 

To derive the total-flux spectra, we extracted all one¬ 
dimensional sky-subtracted spectra optimally (Horne 1986) 
in the usual manner. Each spectrum was then wavelength 
and flux calibrated, and was corrected for continuum atmo¬ 
spheric extinction and telluric absorption bands (e.g., Math- 
eson et al. 2000). The total-flux spectrum was calibrated by 
observing the flux standard star BD-|-17°4708; note, how¬ 
ever, that neither SN 2011dh nor the flux standard were 
observed at the parallactic angle (Filippenko 1982), and so 
the relative-flux calibration of the spectrum is somewhat 
suspect. 

The polarization-angle offset between the half-wave 
plate and the sky coordinate system was determined by ob¬ 
serving the polarized standard star HD 204827 and setting 
the i?-band polarization position angle (i.e., the debi- 
ased, flux-weighted average of the polarization angle over the 
wavelength range 6100-8100 A; see Leonard et al. 2001) to 


© 2013 RAS, MNRAS 000, [T]-?? 







Spectropolarimetry of SN 2011dh 5 


be 59.1°, the value cataloged by Schmidt et al. (1992a). With 
the zeropoint now defined, we also measured the polariza¬ 
tion position angle of the polarized standard star HD 155528, 
and found its value to be within 0.23° of that cataloged by 
Schmidt et al. (1992a). To check for instrumental polariza¬ 
tion, we observed the null standard BD +32° 3739, and found 
it to be null to within 0.07%. 

All data obtained from the Lick, Arizona, and Palomar 
observatories were extracted and calibrated using generic 
IRAI0 routines and our own idl functions. Spectropolari- 
metric analysis was also performed in iraf and idl follow¬ 
ing the methods outlined by Miller, Robinson, & Goodrich 
(1988) and detailed by Leonard et al. (2001) and Leonard & 
Filippenko (2001). All data were corrected for the redshift 
of M51 {z = 0.00155). 


3 ANALYSIS & RESULTS 

Linear polarization is expressed as the quadratic sum of the 
Q and U Stokes parameters, P = ^and the po¬ 
sition angle on the sky is given by 0 = (1/2) tan“^(u/g), 
while carefully taking into account the quadrants in the q-u 
plane where the inverse tangent angle is located. Since P is 
a positive-definite quantity, it is significantly overestimated 
in situations where the signal-to-noise ratio (S/N) is low. 
It is typical to express the “debiased” (or “bias-corrected”) 
form of P as Pdb = {q^ + u^) — (u^ + crj), where ag and 

cr„ are the uncertainties in the q and u Stokes parameters. If 
(cTq + crj) > {q^ + M^), then we set a la upper limit on P of 
In cases where P/ap < 1.5, 9 is essentially un¬ 
determined. All polarized spectra presented herein are dis¬ 
played in this manner. Note, however, that at low S/N, Pdb 
is also not a reliable function, as it has a peculiar proba¬ 
bility distribution (Miller et al. 1988). Thus, for extracting 
statistically reliable values of polarization within a particu¬ 
lar waveband, we have binned the calibrated q and u Stokes 
spectra separately over the wavelength range of interest be¬ 
fore calculating P and 6. All quoted and tabulated values in 
this paper were determined in this manner, although spec¬ 
tra are displayed as Pdb, so they may exhibit slight offsets 
from our quoted values. 

3.1 Weak interstellar polarization 

As light traverses interstellar media, nonspherical dust 
grains aligned with the local magnetic field can polarize inci¬ 
dent photons (Hiltner 1949; Davis 1955). Correcting for this 
component of interstellar polarization (ISP) is often a very 
problematic element of polarimetric analysis, as an accurate 
probe of the ISP from both the Milky Way and, in partic¬ 
ular, the SN host galaxy, can be very difficult to obtain. 
According to Serkowski et al. (1975), the colour excess from 
interstellar extinction can be used to derive an upper limit 
on Galactic ISP, given by Pisp < 9E{B — V)%. The non¬ 
detection of NaD lines at the redshift of M51 by Arcavi et 

^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under a cooperative agreement 
with the National Science Foundation (NSF). 


al. (2011) and Vinko et al. (2012) indicates that extinction 
from the host galaxy is low, with an estimated upper limit 
of E{B — V) < 0.05mag (Arcavi et al. 2011). The Galac¬ 
tic component of extinction in the direction of SN2011dh 
is E{B — V) = 0.035mag (Schlegel et al. 1998), for which 
we expect a maximum possible ISP of < 0.3%, based on the 
upper limit from Serkowski et al. 

To directly gauge the effect of Galactic ISP, we ob¬ 
served the A5 V star HD 117815, which lies < 1° away from 
SN2011dh on the sky and is presumed to be intrinsically 
unpolarized. It has a measured parallax that indicates a dis¬ 
tance of 122 pc (Perryman et al. 1997). At a Galactic lati¬ 
tude of +68°, this corresponds to a Galactic scale height of 
113 pc, so the light from this star should probe the majority 
of intervening dusty ISM from the Galaxy. We observed this 
star with the Kast instrument on two occasions and detected 
no significant polarization, measuring only upper limits of 
Pv < 0.06% on 2011 June 30 and Pv < 0.03% on 2011 July 
6. We also observed the F5V star HIP 65768 with Kast on 
2016 Jun 30. The parallax of this object indicates a distance 
of 211 pc and a Galactic scale height of 154 pc. We measure 
qv = 0.04 ± 0.02% and uv = 0.09 ± 0.02%, which indicates 
P = 0.10+0.03%. At P < 0.1%, it appears that the Galactic 
component of ISP is practically insignificant in the direction 
of SN 2011dh, and much lower than the maximum possible 
value of ISP indicated by the extinction. 

To estimate the vector sum of ISP values from M51 
and the Galaxy, one might proceed by assuming that the 
emission components of the strongest features, such as Ha, 
are completely depolarized, intrinsically (as was assumed in 
the case of SN 1993J; Tran et al. 1997). However, Maund 
et al. (2007a) demonstrated that using the polarization of 
Ha emission for this purpose is probably not reliable, since 
there is the possibility for subtle blending of other features 
(e.g., He I A6678). Furthermore, if a SN has more than one 
component of polarization (for example, Thomson scatter¬ 
ing and/or dust scattering by GSM), then this could pro¬ 
duce net polarization at the same wavelengths as the strong 
intrinsically unpolarized emission features (e.g., see Mauer- 
han et al. 2014). The possible existence of extended GSM 
has been suggested in the case of SN 201 Idh, based on de¬ 
tection of broad Ha in late-time nebular spectroscopy at 
> 200 days (Sahu et al. 2013; Shivvers et al. 2013) and ra¬ 
dio observations (Horesh et al. 2013). Whether this GSM 
component could be influencing the polarization properties 
of SN 201 Idh at the earlier times we are dealing with is 
not clear, but perhaps unlikely, given the lack of substan¬ 
tial UV excess reported for this object (e.g., see Ben-Ami et 
al. 2015). Nonetheless, the possibility of scattering by GSM 
warrants caution in assuming intrinsically depolarized emis¬ 
sion features. Using the Call near-infrared (IR) triplet to 
estimate the total ISP might be a safer approach, since this 
region is devoid of other spectral features. The behaviour of 
this line makes a strong case for very low ISP, since at every 
epoch, q, u, and P all trend close to null (< 0.1%) near the 
wavelength of the emission peak. 

Another method for constraining the total vector sum 
of ISP from the host and the Galaxy is using late-time mea¬ 
surements of the Stokes parameters at the wavelengths of 
any strong forbidden lines as the SN enters the nebular 
phase. Such lines form in low-density regions of the explo¬ 
sion, where the effect of electron scattering is expected to be 
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Figure 1. Temporal evolution of the q and u Stokes parameters (black/blue curves) and total-flux spectra (red curves) of SN 2011dh. 
The data are binned to 25 A, except days 58 and later, which are binned to 50 A. Prominent atomic features are labeled at different 
epochs. Telluric absorption lines in the day-14 data are marked with Earth symbols (circled crosses). 


weak and the polarization should be insignificant. As will be 
demonstrated in the next section, the nebular-phase spectra 
indeed indicate a low value of ISP. We conclude that the 
total vector sum of ISP from the host and the Galaxy in the 
direction of SN 2011dh is not substantial enough to make a 
significant impact on our measurements or scientific inter¬ 
pretation, and we thus move onward, making no attempt to 
remove ISP from the data. 


3.2 Spectropolarimetric evolution of SN 2011dh 

Figure[T] illustrates the evolution of the q and u Stokes pa¬ 
rameters of SN2011dh between 9 and 86 days post explo¬ 
sion, and Figure[2] shows the same evolution for debiased 
polarization P and position angle 6. Each panel in both fig¬ 


ures is accompanied by the total-flux spectrum at the given 
epoch. Table! lists the integrated P-band and continuum 
values of P and 6 for each epoch. The most reliable sample 
of the SN continuum appears to lie in the wavelength range 
7000-7450 A. 

The earliest (day 9) data were obtained when the SN 
photosphere was still in the H envelope, as indicated by 
the Balmer-dominated spectrum. At this stage, the spec¬ 
trum consists of broad P-Cygni lines from Ha, a blend 
of Fe II emission centered near 5300 A with multiple emis¬ 
sion/absorption components extending blueward of 5100 A, 
and a strong feature from the Ca ll IR triplet near 8650 A. A 
relatively weak feature of Na l or He l is also seen in the range 
5800-5900 A, which for a SNIIb is most likely attributable 
to He I A5876 (see Hachinger et al. 2012), and thus we as- 
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Figure 2. Same as Figure[T] but for debiased polarization, P, and position angle, 9. Wavelengths for which Pjcjp < 1.5 have been 
omitted from the 9 curve. Telluric absorption lines in the day-14 data are marked with Earth symbols (circled crosses). 


sume this identification hereafter. At this stage substantial 
polarization is detected. As is commonly seen in CC SNe, 
troughs and crests in P are coincident with the respective 
emission and absorption components of the P-Cygni profiles 
(“inverted P-Cygni” polarization profiles), likely resulting 
from asymmetric or patchy obscuration of the photosphere. 
Regions of broad polarization minima are associated with 
the broad emission components in the total-flux spectra, 
whereby the strong intrinsically unpolarized line emission 
dilutes the polarized flux (Cropper et al. 1988). Fluctua¬ 
tions of ~ 0.5% in P are evident around Hel and Ha, and 
reach ~ 1% for Call IR. A region of line-free continuum 
can be seen in the range 7000-7450 A, where we measure 
P = 0.47 ± 0.03% and 6 = 23.5 ± 1.5°; continuum polar¬ 
ization also probably contributes to the broad hump in P 
centered near 5000 A on days 9 and 14, as indicated by the 


level of P at the continuum near 6050 A. The highest value 
of P is associated with the absorption components of the 
Call IR triplet, reaching ~ 1.1%. Substantial rotations in 6 
are seen across the spectrum, changing by ~ 30° across the 
Fell blend, and fluctuating by ~ 20° across Ha. The largest 
rotations are 80-90° across Call IR and ~ 70° at the far 
blue end of the spectrum, possibly associated with H/9. 

By day 14, Hel features are beginning to appear in the 
flux spectrum. The A5876 line exhibits deepening absorp¬ 
tion, and the A6678 line can be seen as a notch of absorp¬ 
tion near the peak of the Ha emission component. Spec¬ 
tral analysis by Marion et al. (2014) showed that the first 
definitive detection of Hel lines occurred 11 days post ex¬ 
plosion, so at day 14 the photosphere must be entering the 
He layer of the envelope. At this stage, the overall charac¬ 
teristics of P are similar to those on day 9, although the 
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Figure 3. Evolution in the q-u plane of the entire spectrum on days 9, 14, 30, and 36. The data have been binned to 50 A. Colours are 
used to indicate different parts of the spectrum. Note the smaller wavelength coverage on day 14. 


broad enhancement near 6000 A has shifted slightly toward 
the red, with the addition of a narrow peak of enhanced P 
appearing near 6200 A. Substantial evolution in 6 has oc¬ 
curred between days 9 and 14, as the major rotations seen 
across the spectrum have grown in strength, most notably 
near the blueshifted portion Her (on the lower-velocity slope 
of the absorption trough), and also on the red shoulder of 
the emission component. The red end of our wavelength cov¬ 
erage at this epoch cuts off at 7800 A, so the evolution of 
Call IR cannot be ascertained. It should be noted, however, 
that the large rotation in 9 between days 9 and 14 appears 
to be dominated by the evolution of the continuum compo¬ 
nent. The 7000-7450 A continuum region on day 14 exhibits 
P = 0.45±0.02% and 6 = 2.2±1.4°, comparable in strength 
to that on day 9 but shifted by A9 « 20°. 

By day 30, the character of the total-flux spectrum has 
changed considerably. Hel lines have become dominant, in¬ 
dicating that the photosphere has moved into the He layer of 
the outflow. The emission profiles of Ha and Hel A5876 ap¬ 
pear to be asymmetric, having a bump or knee between the 
peak wavelength and the red side of the profile. Enhanced 
P is very prominent for the absorption component of Ha, 
and it exhibits a peculiar fluctuation in strength between 
the blue and red walls of the absorption trough. Specifically, 


P reaches a maximum of ~ 1.2% on the red wall of the 
absorption trough but rapidly drops to null on the bluer, 
higher-velocity side of the trough. Curiously, the associated 
change in 9 across this feature is unremarkable. P drops to 
almost null values near the peaks of He l and Ca ll IR emis¬ 
sion. Less-extreme enhancements in P are also associated 
with the red walls of the absorption components for the var¬ 
ious Hel lines throughout the spectrum. A new feature of 
enhanced P is also seen for the absorption component of O l 
near 7600-7700 A. 

The most prominent rotation is a A9 « 90° shift as¬ 
sociated with Call IR. There is also a ~ 70° rotation as¬ 
sociated with a feature near 7000-7100 A, which might be 
attributable to Hel A7065, although in this case it would 
be unclear why similar rotations are not also seen for other 
Hel features. The continuum region used for days 9 and 14 
appears to have dropped to P = 0.18 ± 0.04% and rotated 
to 9 = 12.0° ± 3.4°, but this could be heavily influenced by 
the new appearance of line structure. On the other hand, 
the broad polarization hump centered near 5000 A on days 
9-14, which is continuum polarization, has dropped as well. 
It thus appears that the continuum polarization decrease 
between days 14 and 30 could be real and substantial. 

The day-36 spectrum exhibits only slight changes in 
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Figure 4. Evolution in the q-u plane for prominent spectral lines, colour coded by radial velocity. Four epochs are shown for Ho and 
He I A5876, and three epochs are shown for the Call IR triplet and Ol. The inset frames display the flux spectra with arbitrary scaling. 
The dashed line represents the “axis” defined by the day-9 continuum (black triangle). The Stokes data are binned to 50 A. The black 
error bar in each panel represents the average uncertainty. 


Ha and Hel polarization with respect to day 30, and the 
emission profiles still appear to exhibit the asymmetric knee 
on the red sides. However, substantial changes in P and 
9 are evident for other lines. In particular, we see strong 
enhancement in P for the highest-velocity absorption of Ca ll 


IR. Enhanced P associated with Hel A5876 and O l have also 
grown to 1%. 

By day 65, the absorption component of Ha is almost 
completely gone, implying that the H-rich outer layers of 
the outflow are close to becoming fully transparent. P-Cygni 
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profiles of Hel A5876, Can IR, and Ol are still prominent, 
however, each maintaining enhanced P levels of ~ 1% for 
their blueshifted absorption components. Forbidden emis¬ 
sion lines of [Fell] and [Call] have appeared in the wave¬ 
length range 7000-7400 A, heralding the nebular transition. 
By day 86, the P-Cygni absorption components of all lines 
have weakened as the explosion thins out, but strongly en¬ 
hanced P is still seen for blueshifted Hel A5876 and for Can 
IR, in addition to weaker yet significant features across the 
rest of the spectrum. We note that the asymmetric knee on 
the red side of the Hel A5876 profile, which first became 
apparent on day 30, is still present through day 86. 

Since the forbidden lines seen at late times form in low 
density regions of the outflow, electron scattering is weak, 
and thus the polarization should be insignificant near the 
central wavelength of each line. As an additional means of 
estimating the ISP, we averaged the Stokes parameters over 
the regions corresponding to the peak fluxes of the lines 
(7020-7110 A and 7270-7350 A). We obtained qisp = 0.08 ± 
0.05% and uisp = —0.01 ± 0.05%, which implies that P = 
0.08 ± 0.07% at these wavelengths. The low value of ISP 
indicated by this result supports our decision in §3.1 to make 
no attempt at removing ISP from the data. 

3.2.1 Behaviour in the q-u plane 

Figure^] shows the q-u plane evolution over the full wave¬ 
length range of the spectra for days 9, 14, 30, and 36. Fig- 
ure[3]illustrates the behaviour in the q-u plane across promi¬ 
nent line features, and also shows the flux spectra of the line 
plotted with respect to radial velocity (RV). On day 9, the 
full range of data forms a cloud of points that are roughly 
centered in quadrant:. The large excursion near the Can 
IR triplet is clearly seen as a looping pattern that moves far 
into quadrant ll, with maximum excursion occurring for RVs 
ranging from —10, OOOkms”^ to —15, 000kms“^. This fea¬ 
ture is mostly nonaxisymmetric relative to the 7000-7450 A 
continuum sample (shown as a dashed line). However, the 
fastest Can absorption on day 9 has RV = —20, 000 kms~^ 
and moves back into quadrant l, nearer to the continuum 
axis. 

Blueshifted Ha absorption creates a loop in the q-u 
plane that, on day 9, deviates from the continuum axis for 
RV below —20, 000 kms“^. On day 14, the large rotation of 
9 near Ha is clearly apparent in the q-u plane, but by days 
30-36 the q-u loops associated with this line, and with Hel, 
become more elongated and roughly aligned with the posi¬ 
tion angle of the day 9-14 continuum axes. The velocities at 
maximum polarization for H and He are u ~ —10, 000 km s“^ 
and V ~ —8000kms“^, respectively. The pronounced varia¬ 
tion in 9 observed across Ha on day 14 could be the result 
of the evolving continuum polarization level; the changing 
vector addition of the Stokes parameters of this continuum 
component with those of Ha could produce an artificial ro¬ 
tation across the line, rather than being, for example, an 
intrinsic change in the geometry of Ha absorption itself. 

The Ol feature, like Can IR, loops in a highly nonax¬ 
isymmetric manner with respect to Ha, Hel, and the day-9 
continuum (most apparent on days 30-36), although the q-u 
loop of this line drifts into quadrant iv rather than quadrant 
II, as Call IR does. The velocity at maximum polarization 
for Call is comparable to that of Ha, while for O l the veloc- 
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Figure 5. Rotated stokes parameters for SN 2011dh on days 9 
and 30, aligned to the axis of symmetry defined by the day-9 
continuum region (7000-7450 A). Total-flux spectra are shown in 
red. 


ity is —4000kms“^, substantially slower than Call, Ha, 

and Hel. 

To further illustrate the distribution of line features 
with respect to the day-9 continuum axis, we rotated the 
q and u Stokes parameters (RSP) to 9 = 23.5°; the results 
are shown in Figure[5l It can be seen that the majority of 
polarization at wavelengths shortward of 8000 A lies in qnsp, 
along the axis of the continuum. The q-u loop created by 
the Call feature is evident as a large deviation in the RSP, 
with most of the polarization contained in mrsp (nonaxisym¬ 
metric with the continuum). We examine this feature more 
closely in Figure[6l which illustrates qnsp, mrsp, and flux 
as a function of RV for Call and Ha. It can be seen that 
on day 9, the polarized profile of Call is relatively broad 
compared to later epochs. In the flux spectrum there ap¬ 
pears to be two distinct components of absorption, with the 
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Figure 6. Evolution of qrsp (thin black curves) and tiaSP (thick 
black curves) as a function of radial velocity for the Ca 11 IR triplet 
and Ha. Total-flux spectra are shown in red. In the upper-right 
panel, the day-14 data are represented by the dashed curves. 


most highly blueshifted component exhibiting an edge that 
extends to RV « —20, 000 kms”*^, and the same is true for 
Ha. For Call, most of the polarized flux on day 9 is con¬ 
tained in mrsp for RV < 15,000 kms“^, while at higher RV 
some of the polarization shifts into gasp. By day 36, the 
polarized feature exhibits a lower range of velocities cen¬ 
tered at ~ —12, 000 kms“^. Note that the emission features 
of Ca II and Ha both exhibit a significant blueshift of 2000- 
3000 kms“^ on day 9, and appear to become more centered 
near rest velocity by days 14-30. 


4 DISCUSSION 

4.1 Sources of continuum and line polarization 

4-1.1 The potential effect of stellar deformity 

The substantial continuum polarization detected from 
SN 2011dh at the earliest epoch on day 9 implies that the 
photosphere exhibited global aspherlcity while in the outer 
H-rich envelope. When interpreted in terms of the oblate, 
electron-scattering model atmospheres of Hoflich (1991), our 
results imply an asymmetry of at least 10%; the actual value 
naturally depends on the (unknown) viewing orientation, 
and on the extension and depth of the electron-scattering 
atmosphere. The subsequent drop in continuum polariza¬ 
tion to P « 0.2% by day 30, if real and not the result of 
line contamination starting at that epoch, suggests that the 
photosphere took on a more spherical configuration after 
passing into the He and C-l-0 layers. This behaviour is op¬ 
posite the polarimetric trends observed in some well-studied 
SNell-P, which begin with low values of P during the first 
few months, indicative of a mostly spherical outer envelope, 
and then increase in polarization only after H recombines 
and the photosphere enters the He and C-l-O layers; in these 


cases, the global aspherlcity has its origin in the core (Jeffery 
1991; Leonard et al. 2001, 2006; Wang & Wheeler 2008). 

For SN 2011dh, a higher degree of global aspherlcity 
on the outside possibly suggests that some other physical 
mechanism, perhaps related to binary evolution, might be at 
work. After all, in the majority of cases, SNIIb progenitors 
are thought to have lost the bulk of their H envelopes to tidal 
stripping by a massive accretor (Podsiadlowski et al. 1993; 
Hoflich 1995; Smith et al. 2011; Dessart et al. 2012). Inter¬ 
estingly, a hot star has now been detected at the explosion 
site of SN 2011dh, and it is a viable candidate companion to 
the progenitor (Folatelli et al. 2014). If the strong gravita¬ 
tional tides that induce mass transfer can also geometrically 
distort the envelope of the progenitor star, or synchronise 
the star’s rotation at a fast enough rate to oblately deform 
it, then the star could develop an aspherical density profile. 

Interestingly, recent interferometric observations of the 
interacting yellow hypergiant binary HR5171A have pro¬ 
vided evidence for significant envelope distortion in the pri¬ 
mary (Chesneau et al. 2014). Although more luminous than 
the progenitor of SN 2011dh, the example of HR5171A in¬ 
deed motivates the expectation of binary-induced stellar dis¬ 
tortion. If similarly induced geometric deformity in the pro¬ 
genitor of SN 2011dh is substantial enough, the resulting as¬ 
pherical density gradient could provide a means to generate 
an aspherical shock front and explosion geometry (Yamada 
& Sato 1990; Steinmetz & Hoflich 1992). It should first be 
explored, however, whether sufficient tidal deformation is 
plausible for the SN 2011dh progenitor system. 

To assess the potential for tidal deformation, we require 
estimates of the physical parameters of the binary system 
and its component stars. For the case of SN 2011dh, Ben¬ 
venuto et al. (2013) performed binary stellar evolutionary 
calculations that followed the simultaneous evolution of both 
stars in the system to produce a pre-SN progenitor that is 
consistent with the space-based photometry (Van Dyk et al. 
2011, 2013). Benvenuto et al. concluded that the most plau¬ 
sible model is one in which the primary and secondary com¬ 
ponents of the binary began with respective initial masses 
of Ml = 16 Mq and = 10 Mq, and an orbital period of 
P — 125 d. It had a moderate mass-transfer efficiency fac¬ 
tor oi P — 0.5 and a specific angular momentum value of 
a = 1 (i.e., the angular momentum of the material lost by 
the system is equivalent to the angular momentum of the pri¬ 
mary star). At the moment of explosion, the model predicted 
that the YSG primary would have a heavily reduced final 
mass of M/ « 4 M©, a small hydrogen content of Mh ~ 3- 
5 X 1O“®M0, and a final radius of Rysg ~ 250 R© (con¬ 
sistent with the photometry); the mass-gaining secondary 
ended up with a final mass of M/ « 16 M©. Note that /? is 
one of most unconstrained parameters in binary evolution; 
regardless, considering the full range of possible values in 
their models (0 < /3 < 1), Benvenuto et al. concluded that 
the final semimajor axis of the system at the time of ex¬ 
plosion lies in the range Af = 842-1009 R©, which implies 
A//i?YSG ~ 4. Meanwhile, the ratio of the YSG’s Roche- 
lobe radius and the final separation can be approximated by 
Rnc/Af « 0.46Mi/(Mi-|-M 2)^^^ (Eggleton 1983), which in 
our case yields Rrl/M/ « 3.5. Thus, Rysg ~ Rrl. 

Frankowski & Tylenda (2001) present analytic solutions 
for the geometric dimensions of a triaxial ellipsoid (TE) that 
stars will develop in response to tides in a massive binary (a 
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football-like shape). The principal axes of the ellipsoid are 
defined by the direction toward the companion (semi-axis 
ai), the direction in the orbital plane that is perpendicular 
to the latter (semi-axis 02 ), and the direction perpendicular 
to the orbital plane (semi-axis a^), each given by 


ai = R 


1 - 1 - 





a2 = R 


1-f 



(1 



as = R 1 



(2 + 5q) 



( 1 ) 


where q is the mass ratio, R/A is the ratio of stellar radius 
to binary separation, and 7 is the “gravity effectiveness” 
parameter, which is some value between zero and unity. For 
relatively compact main-sequence dwarfs, 7 « 1 is an ac¬ 
ceptable approximation that leads to only very subtle de¬ 
formity for R/A « 1/4. For YSGs, which are dynamically 
unstable (e.g., see Stothers & Chin 2001), 7 could be ~ 0.1 
or lower, as a result of the relatively intense radiation pres¬ 
sure on the envelope, strong winds and high mass-loss rates, 
convection, and pulsations; such factors will profoundly in¬ 
fluence the effective gravity and Roche potential, and can act 
to shrink the effective Rrl. If we therefore assume 7 = 0.1 
and adopt the values of q and R/A above. Equation (9) 
from Frankowski & Tylenda (2001) implies TE semi-axes of 
(ai : 02 : 03 ) ~ (1.1: 1.0 : 0.9), or ~ 10% deviation from spher¬ 
ical symmetry along two axes. Lower values of 7, q, and R 
will, of course, result in a higher degree of deformity, as will 
smaller values of A. 

Stellar asphericity can also result from rotation, which 
might become fast enough to deform the star into an oblate 
geometry if the binary is tidally locked; this is a likely pos¬ 
sibility for extended stars that nearly fill their Roche lobes. 
If so, the binary parameters adopted from Benvenuto et al. 
(2013) imply that the YSG would have a rotational veloc¬ 
ity of vysg ~ 17kms“^. The critical rotational velocity, at 
which point the star will begin to break up, is given by 
Ucrit = 7 /GMett/R (Langer et al. 1997; Langer 1998), where 
Meff = M(1 - Fe) and F^ « 2.1 x 10"® L/M. For the YSG 
progenitor, the values L = lO"^'® L© and M/ = 4 M 0 im¬ 
ply Meft = 3/4M/ = 3 M 0 ; thus, Ucrit = 48kms“^ and 
VYsc/'i^crit ~ 0.35. This subcritical value would generate a 
rotationally induced asphericity of < 5%. Of course, this 
all depends on the model-dependent binary parameters we 
have adopted; a shorter binary separation or larger stellar 
masses would result in faster stellar rotation for a tidally 
locked system, and thereby larger asphericities. 

A deformed stellar envelope, whether induced by tidal 
stretching, rotation, or both, will exhibit an aspherical den¬ 
sity profile, having a steeper gradient along the poles; this 
can have a profound effect on the geometry of the pho¬ 
tosphere of the subsequent SN explosion. Specifically, hy¬ 
drodynamic simulations of shock propagation in rotating 
blue supergiant stars (like the progenitor of SN 1987A) have 
demonstrated that an initially spherical shock wave that 
encounters a steeper density gradient along the poles will 
develop a prolate-shaped front and break out of the poles 


first (Steinmetz & Hoflich 1992). The interaction of this pro¬ 
late shock with the oblate density stratification in the star 
causes a second shock wave that moves from the poles to the 
equator. The resulting transport of energy and momentum 
subsequently produces a quasihomologous outflow having an 
oblate outer and prolate inner density stratification. In their 
simulations, Steinmetz & Hoflich tracked the evolution of 
the observed polarization from the electron-scattering pho¬ 
tosphere, predicting a relatively high continuum polarization 
at early times, while the photosphere is in the outer enve¬ 
lope {P « 0.4-1.5% during the first 10 days), followed by a 
rotation in 6 and diminishing P as the photosphere passes 
into the respective prolate and spherically stratihed layers 
of the interior. 

Although the detailed outcome depends on the stellar 
radius, detailed density profile, rotation (differential versus 
rigid), and inclination angle, the basic predictions of the 
above simulations appear to be more or less consistent with 
the polarized continuum evolution of SN 2011dh, and thus 
might provide some general insight. The future construc¬ 
tion of similar hydrodynamic models of shock propagation 
in rotating and/or deformed YSGs could help put the hy¬ 
pothesis of a connection between progenitor deformity and 
SNIIb explosion geometry on a more quantitative footing, 
leading to more specific predictions for how particular spec- 
tropolarimetric observables and their associated timescales 
can constrain the physical properties of SNellb and their 
progenitors. 


4 . 1.2 Oblique shock breakout 

Recent calculations by Matzner et al. (2013) have demon¬ 
strated how nonradial fluid motions generated by aspheri¬ 
cal core collapse can dramatically alter the dynamics and 
emission of the emergent shock. Depending on the degree 
of asphericity, the shock front can develop an oblique orien¬ 
tation with respect to the stellar surface over a signihcant 
fraction of the star. The subsequent outflow from regions 
crossed and energised by this shock becomes unsteady, po¬ 
tentially stifling the intensity and duration of the breakout 
flash, limiting ejecta speeds, and casting matter sideways. 
Gollisions of matter outside the star might also occur and 
significantly influence the observational features of the re¬ 
sulting transient. 

Matzner et al. (2013) discussed SN2011dh in the con¬ 
text of oblique shock breakout, and suggested that the re¬ 
sulting asphericity of the highest-velocity ejecta might ex¬ 
plain why the early-time luminosity and radio emission indi¬ 
cated a compact progenitor (Arcavi et al. 2011; Soderberg et 
al. 2012 ), even though post-explosion images clearly showed 
that the SN stemmed from an extended YSG (Van Dyk et 
al. 2013); the reason for the discrepancy lies in the unrealis¬ 
tic assumption of spherical symmetry. Interestingly, oblique 
shock breakout in SN 2011dh might also explain why the SN 
exhibited the strongest continuum polarization during the 
earliest stages, since the photosphere was in the outer enve¬ 
lope of the quasihomologous flow where velocities are high¬ 
est. This possibility motivates the construction of hydrody¬ 
namic models that will yield information on the aspherical 
explosion geometries that stem from oblique shock break¬ 
out, enabling the further development of radiative transfer 
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and scattering models that will predict the polarization evo¬ 
lution for such a scenario. 

4-. 1.3 Asymmetric radioactive heating 

Net continuum polarization from a nonax;isymmetric pho- 
tospheric configuration could also result from an aspheri- 
cal distribution of radioactive heating and ionization. For 
many other CC SNe, clumps or jet-like configurations of 
®®Ni rising rapidly from the core have been invoked to ex¬ 
plain their spectropolarimetric properties at various phases 
(Chugai 1992; Leonard et al. 2001, 2006; Chugai et al. 2005; 
Maund et al. 2007a, 2007b; Chornock et al. 2011), and their 
asymmetric nebular line proHles (Filippenko et al. 1994; Spy- 
romilio 1994; Matheson et al. 2000; Milisavljevic et al. 2010; 
Chornock et al. 2011; Smith et al. 2015), including the late¬ 
time nebular spectrum of SN2011dh (Shivvers et al. 2013). 

The hypothesis of aspherical radioactive heating is well 
motivated; indeed, modern three-dimensional (3D) hydrody¬ 
namic models of CC explosions have shown that relatively 
small convective asymmetries initiated by fast-rising plumes 
of neutrino-heated matter in the core can, within hours, 
grow into large-scale asymmetries in the heavy-element dis¬ 
tribution of the outflow (Hammer et al. 2010; Wongwatha- 
narat et al. 2010). Models for both red and blue supergiant 
progenitors with 10 M© H envelopes produce ®®Ni plumes 
that mix out to velocities up to 3000-5000 kms“^ in the ho¬ 
mologous flow (Hammer et al. 2010). More relevantly, 2D 
hydrodynamic models of SNellb with H envelopes that are 
several orders of magnitude less massive indicate that bub¬ 
bles or plumes of ®®Ni could mix farther out to 6500 km s“^ 
(Iwamoto et al. 1997), near the H/He interface (one might 
expect the mixing to extend to even larger radii in future 3D 
models of SNellb, if the relative velocity predictions from 
the aforementioned 2D and 3D models of H-rich CC SNe are 
mirrored in SNellb). 

Rising ®®Ni plumes can snowplow surrounding material 
into a shell that can become Rayleigh-Taylor (RT) unstable, 
potentially mixing ambient envelope material into the bub¬ 
bles as they rise (Basko 1994). Depending on the path length 
for the thermalization of gamma rays from the decay of ®®Ni 
or ®®Co, radioactive plumes will be encased in cocoons of en¬ 
hanced free-electron density and excess Thomson scattering 
(e.g., see Chugai et al. 2005). If the cocoons extend into the 
H layer, they could produce an irregular or mottled pho- 
tospheric geometry at early times, thus contributing to the 
net continuum polarization (e.g., on days 9 and 14). Depend¬ 
ing the projected position of a given plume relative to the 
global underlying photosphere, the local polarization inte¬ 
grated over the plumes spheroidal surface could effectively 
cancel out, having essentially the same effect as a photo¬ 
sphere that is partially obscured by line absorption. This 
interpretation potentially explains the rough alignment be¬ 
tween the position angle of the early day-9 continuum axis 
and that of H and He l line polarization on days 30-36, which 
is further illustrated by Figure[71 In this scenario, we sug¬ 
gest that the same spheroidal plumes that generate locally 
enhanced Thomson scattering at early times (and where 
the integrated polarization locally cancels out because of 
the spherical pattern on the plume) later become plumes of 
relatively high line optical depth for Ha and Hel absorp¬ 
tion, since the enhanced electron density in those regions 
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Figure 7. Evolution of P and 9 for Ha from days 30 and 36. 
The value of 9 for the day-9 continuum is represented by the 
horizontal grey line in the lower panel. Spectral regions of interest 
are highlight with bold lines. 


also helps to populate the upper atomic levels of the transi¬ 
tions (e.g., Lucy 1991; Fassia et al. 1998; Chugai 1992, 2006; 
Chugai et al. 2005; Swartz et al. 1993; Hoflich, Khokhlov, & 
Wang 2001). Thus, the plumes at early times will break the 
symmetry of the global photosphere in a geometric configu¬ 
ration that is similar to symmetry breaking by line absorp¬ 
tion at later times, in both cases producing net polarization 
with similar position angles. Of course, the construction of 
detailed 3D scattering models that account for the geome¬ 
tries discussed above are needed to put this hypothesis on a 
more quantitative footing. 

An aspherical distribution of plumes in the ejecta might 
also explain the knee feature on the red sides of the Hel 
A5876 and Ha emission profiles observed on day 30 and later, 
which could be an imprint of line emission from such struc¬ 
tures on the receding hemisphere of the explosion, but just 
outside the limb of the projected photosphere on the sky. 
Interestingly, similar morphology was seen in the Ha line of 
SN 1987A on day 71, interpreted as the imprint of clumpy 
ejecta on the receding hemisphere of the explosion (Chugai 
1992). In that case, moreover, Thomson scattering of pho- 
tospheric radiation by radioactive plumes on the receding 
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Figure 8. Illustration of polarized line formation in SN 2011dh (not to scale). The formation of the polarized Ha profile is represented 
here (inset at lower left, also showing flux as thinner lines), but the basic idea is applicable to other chemical species that exhibit enhanced 
line polarization. The upper diagrams are face-on views, and the lower diagram illustrates the transverse view. In the homologous flow, 
isovelocity planes (dashed lines) are perpendicular to the line of sight. In the face-on views above, isovelocity planes having the largest 
radial velocity slice through a smaller “cap” of the sphere closest to the observer, which is why the circular region encompassed by 
the dashed line is smaller for higher blueshift. The fastest —16, 000 km s“^ isovelocity plane (green dashed line) intersects the patchy 
tenuous outer regions of the flow, which asymmetrically obscures the electron-scattering photosphere (arrow-marked region) and give 
rise to net polarization at those radial velocities. A relatively smooth or tenuous region of the hydrogen envelope is intersected by 
the —14, 000 kms“^ isovelocity plane (blue dashed line), resulting in uniform or weak obscuration of the photosphere and, hence, low 
polarization. The —10, 500kms“^ isovelocity plane (red dashed line) intersects several high optical depth (t) regions of Ho/He I, where 
“cocoons” of enhanced excitation for these transitions form around plumes of radioactive ®®Ni, which unevenly obscure the underlying 
photosphere, giving rise to the peaks in the polarized line profiles. 

hemisphere of the explosion was successfully modeled as the 
source of polarization in SN 1987A. 

FigureO illustrates a configuration that potentially ex¬ 
plains the morphology of the polarized line profiles in the 
quasihomologous outflow. Homology implies that an imagi¬ 
nary plane intersecting the flow perpendicular to our line of 
sight (in the plane of the sky) will sample a region of con¬ 
stant RV over that entire plane; we refer to these imaginary 
cross sections as “isovelocity planes.” The peak RV of line 
polarization is a useful indication of the depth of patchy ab¬ 
sorption by clumps in the quasihomologous flow, and can be 
used to probe the relative geometries of atomic species at 
various radii (Kasen et al. 2003). For example, the fastest 


outermost isovelocity planes with RV « —16, 000 km s“^ in¬ 
tersect patches of high-velocity ejecta in the outermost re¬ 
gions of the flow, possibly consisting of clumps from the 
progenitor star’s wind swept up and accelerated by the ex¬ 
plosion. These patches partially obscure the polarized pho¬ 
tosphere and produce enhanced polarization at those corre¬ 
sponding RV values. Deeper in the homologous flow, isove¬ 
locity planes with RV « —14, OOOkms”^ intersect a rela¬ 
tively smooth or tenuous region in the thin H envelope that 
obscures the photosphere evenly or weakly and thus pro¬ 
duces no net polarization (creating the steep drop and local 
minimum of the polarized profile near this value of RV). 
Deeper isovelocity planes with RV « —10, 500kms~^ inter- 
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Figure 9. Left-. Enhanced polarization as a function of radial velocity for Ho, He I, Call, and Ol on day 36 (and Ho on day 14; dotted 
line). Homologous expansion is indicated by the relatively fast velocity for Ho and Can, compared to the intermediate velocity for He I 
and the slow velocity for Ol. Right: Temporal evolution of the radial velocity for the polarized line peaks. 


sect the clumpy region produced by the ®®Ni plumes, cre¬ 
ating the peak of the polarized Ha and Hei line profiles at 
this RV. 


4-1-4 Outer Call and inner Ol 

As illustrated by Figure|9l Ha and Call absorption on day 
9 exhibit the highest maximum RV observed in SN2011dh. 
At this time, the emission peaks of both lines have rela¬ 
tively high blueshifted velocities of 2000-3000 km s“^, with 
Ha moving back toward rest velocity by day 14 (note that 
the wavelength range of the day-14 data does not sample 
Call). Such early blueshifts are a common feature of SNell; 
the phenomenon has been attributed to the steep density 
profile of the outflow, which results in highly localised emis¬ 
sion regions and exacerbates obscuration of the receding 
hemisphere of the explosion, resulting in a net blueshift of 
emission peaks during the first month or so (e.g., see Ander¬ 
son et al. 2014, and references therein); the relatively low- 
H envelope masses in SNellb might explain why the early 
blueshift lasts for < 14 days. It thus seems plausible that 
highest-velocity Call in the outermost regions of the flow 
is, like high-velocity Ha, part of the progenitor’s outer en¬ 
velope or wind swept up and accelerated by the explosion. 
The apparent antisymmetry of the polarized Call absorp¬ 
tion feature with respect to the day-9 continuum (see Fig- 
ures[5] and [6ll is consistent with this interpretation, demon¬ 
strating that the patchy absorption on day 9 does not share 
the geometric configuration of the continuum opacity source 
or the prominent polarized features of Ha and He l on later 
days. Meanwhile, the disappearance of the high-velocity ab¬ 
sorption components and the shift of emission peaks back 
toward rest velocity indicates that this separate outer re¬ 
gion has become transparent by day 14. 

O I exhibits the lowest RV and a significantly different 


position angle than the other lines, and thus appears to be 
the hallmark of core oxygen synthesised by the progenitor 
star during post-main-sequence evolution. Interestingly, the 
RV of O I is consistent with the results of modern 3D hydro- 
dynamic simulations of CC explosions, which predict core 
O to be mixed out to the 4000-5000 km s“^ layers of the 
homologous flow (Hammer, Janka, & Muller 2010). More¬ 
over, the different position angle indicates that the inferred 
dumpiness of O l could have a different physical origin than 
that of Ha and Hel, perhaps arising from fluid instabili¬ 
ties deeper in the outflow. For SNellb with low-mass H en¬ 
velopes, large-scale RT instabilities are unable to develop 
at the H/He interface, since there is an insufficient overlying 
mass to facilitate their growth. RT structure is, however, ex¬ 
pected to develop at the deeper He/C-fO interface (Iwamoto 
et al. 1997), and thus might be another potential source of 
patchy OI absorption. 


4-1.5 The potential effect of scattering by CSM 

Another possibility that should be considered for producing 
continuum polarization is scattering of SN photons by CSM. 
If the CSM is sufficiently dense, strong interaction will give 
rise to an additional vector component of Thomson scatter¬ 
ing near the interface of the SN shock and CSM (Hoffman 
et al. 2008; Mauerhan et al 2014); in the case of dense CSM, 
however, one might expect the total-flux spectrum to de¬ 
velop narrow emission-line profiles from CSM that becomes 
photoionized by the UV/X-ray emission at the shock/CSM 
interface (i.e., this should produce a SNIIn spectrum). No 
such narrow lines were seen in our flux spectra of SN 2011dh 
on day 9, or in spectra from others obtained as early as day 
4 (Marion et al. 2014), so Thomson scattering from dense 
ionized CSM seems to be an unlikely alternative for pro¬ 
ducing an aspherical electron-scattering pseudophotosphere 
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during the first month of our spectropolarimetric data. How¬ 
ever, Thomson scattering from CSM interaction cannot be 
ruled for the earliest phases of SN 2011dh, where spectropo¬ 
larimetric data are not available. Indeed, the very early ra¬ 
dio and millimeter detection of SN2011dh just 3 days post 
explosion indicates substantial interaction between the SN 
shock the extended stellar wind of the progenitor (Horesh et 
al. 2013). Recent “flash spectroscopy” of other CC SNe, from 
just several days to as early as hours after explosion, have 
indeed exhibited evidence for short-lived CSM interaction 
and Thomson scattering via Lorentzian emission-line profiles 
(e.g., Gal-Yam et al. 2014; Shivvers et al. 2015; Smith et al. 
2015). Therefore, it is plausible that SN2011dh produced a 
detectable spectropolarimetric signal from CSM interaction 
before our earliest day-9 data. 


Alternatively, small particles of dust in CSM could re¬ 
sult in Rayleigh scattering of the optical SN photons (i.e., a 
“light echo” effect; Chevalier 1986), which will linearly po¬ 
larize a substantial fraction of the reflected light. As argued 
by Wang & Wheeler (1996) in the case of SN 1987A, if the 
dust sublimation timescale is longer than the duration of the 
UV flash from the SN, small surviving grains can produce an 
effective medium for optical scattering. This was considered 
to be a very likely explanation for the strong continuum po¬ 
larization detected from SN 1993J (Hoflich et al. 1996; Tran 
et al. 1997). As mentioned above, early radio and millime¬ 
ter observations of SN 2011dh provided direct evidence for 
CSM from the progenitor wind (Horesh et al. 2013), while 
late-time X-ray emission roughly 500 days after explosion in¬ 
dicated the presence of dense CSM at ~ 8 x 10^® cm (Maeda 
et al. 2014), also consistent with a steady-state YSG wind. 
The corresponding light travel time across this radius is ~ 31 
days. Therefore, continuum polarization via optical scatter¬ 
ing is plausible for the first month of SN2011dh if enough 
dust grains in the wind survived the UV flash of the SN. 
However, in this case the apparent drop in the continuum 
polarization near the spectral transition from H-rich to He- 
rich would have to be a simple coincidence unrelated to the 
photospheric evolution. 


A potential factor that could help distinguish between 
Rayleigh scattering by small dust particles and Thomson 
scattering by free electrons is the wavelength dependence 
of P for the former case, owing to the wavelength depen¬ 
dence of the Rayleigh scattering phase function. For Thom¬ 
son scattering, on the other hand, P is wavelength indepen¬ 
dent. A Rayleigh continuum component, if present, will be 
most prominent at wavelengths < 4500 A. The day 14 data 
from Bok/SPOL provide coverage in the range 4000-4500 A, 
where the average polarization is P — 0.09% + j — 0.03%, 
nearly null and substantially lower than the continuum level 
at this epoch. This appears to be inconsistent with Rayleigh 
scattering by dusty CSM. Instead, the observed character¬ 
istics seem consistent with the “depolarizing” effect of line 
blanketing in this region of the spectrum by material above 
the Thomson-scattering photosphere. We therefore conclude 
that the day 9-14 continuum polarization probably does not 
include substantial contribution from Rayleigh scattering by 
dusty CSM. 


4.2 Comparison to other SNellb 

Despite their relatively low rates among CC explosions, 
SNellb are becoming relatively well represented in terms 
of objects for which there is high-quality spectropolarime- 
try. Although thorough temporal sampling is often not ob¬ 
tained, the growing sample is beginning to facilitate a search 
for systematic spectropolarimetric trends that might point 
to common physical scenarios for SNellb and other classes 
of SNe as well. 

Figure[Tn] compares the total-flux spectra and polariza¬ 
tion of the ellb SNe2011dh (day 30) and 1993J (day 33; 
Tran et al. 1997), and the cllb SNe2001ig (day 31; Maund 
et al. 2007a) and the Ib—J-IIb transitional SN 2008ax (day 
14; Chornock et al. 2011), all near phases where Hel has 
recently become prominent in the spectrum. SN 2011dh and 
SN 1993J exhibit the maxima of their enhanced Hq and He l 
line polarizations at the same wavelengths and, hence, the 
same velocities, while SN 1993J exhibits slightly stronger 
values of polarization for Ha and substantially stronger po¬ 
larization for Hel. Relatively strong continuum polarization 
is also evident for SN 1993J in several regions of the spec¬ 
trum, while the continuum polarization of SN2011dh and 
SN 2008ax has diminished before the phases represented in 
Figure lTOl but has yet to rise in the case of SN 2001ig (see 
Maund et al. 2007a). SN 2008ax stands out as having much 
stronger peak and higher RV Ha polarization than the oth¬ 
ers, by a large margin. 

The conjecture that most SNellb exhibit substan¬ 
tial continuum polarization at early phases (typically near 
0.5%), relative to the frequent absence of continuum polar¬ 
ization in the early stages of SNell-P, is an interesting trend 
that could help to reveal important clues about the nature 
of these explosions and their progenitors. Indeed, if most 
SNellb stem from close interacting binary systems, then the 
imprint of tidal deformation on the continuum polarization 
(see §4.1.1) might be present in many cases. Tidal effects 
were suggested for SN 1993J to potentially explain its in¬ 
ferred ellipsoidal photosphere (Hoflich et al. 1996; Trammell 
et al. 1993). Indeed, the progenitor had a radius that was 
~ 3 times larger than that of SN 2011dh, so the influence of 
tidal distortion on the spectropolarimetry could potentially 
be more extreme and longer lasting in that case. For the cllb 
subset of SNe that stem from relatively compact progenitors 
(perhaps, stripped-envelope Wolf-Rayet stars), direct tidal 
stretching of the progenitor envelope might seem less likely, 
unless the binary separations are relatively small. If cllb pro¬ 
genitors are relatively fast rotators, however, then the result¬ 
ing aspherical density profile for the star could also lead to 
a similar geometric evolution for the SN shock and the sub¬ 
sequent photosphere (direct tidal deformation by a binary 
companion, or fast rotation, will both give rise to a steeper 
density gradient along the polar axis of the star; see §4.1.1). 
Interestingly, the transitional cllb SN 2008ax exhibited sig¬ 
nificant continuum polarization on days 6 and 9 (P « 0.6%), 
indicating that the outer envelope was indeed quite aspher¬ 
ical; rotational deformation thus might be a plausible ex¬ 
planation in this case. For the separate case of the cllb 
SN 2001ig, however, the evolution from relatively weak con¬ 
tinuum polarization during the H-rich phase (P « 0.2%) to 
a strong value of ~ 1% after the photosphere passed into 
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Figure 10. Comparison of the total-flux spectra and RSP (ro¬ 
tated by constant angles deflned by continua axes) for the ellb 
SNe2011dh (day 30; black) and 1993J (day 33; red; Tran et al. 
1997), and the cllb SNe2008ax (day 14; blue; Chornock et al. 
2011) and 2001ig (day 31; green; Maund et al. 2007). 


the He layer (Maund et al. 2007a) is behaviour more in line 
with observations of SNell-P. 

Alternatively, the effects of aspherical radioactive heat¬ 
ing by jets or plumes of ®®Ni might also explain the com¬ 
monality of continuum polarization in SNellb. The earlier 
appearance of this effect, relative to SNell-P with massive 
H envelopes, is naturally explained by the relatively low- 
mass envelopes of SNellb. Furthermore, the apparent con¬ 
nection between the position angles of the early-time contin¬ 
uum, Ha, and Hel line polarizations in SN 2011dh, if real, 
might provide an additional reason to favour this physical 
scenario (see §4.1.3). The same connection was proposed 
for SN 2008ax, but in that case the similarity between the 
line and continuum position angles was even more obvious 
(the system exhibited a higher degree of axisymmetry) and 
was shared by the Ol line as well (Chornock et al. 2011). 
The possibility that SN2011dh and SN 2008ax are similar 
in this regard, even though they are members of the dif¬ 
ferent ellb and cllb subclasses, might provide a reason to 
favour aspherical radioactive heating as the cause of their 
early-time continuum polarization. After all, the origin of 
radioactive plumes is in the core, and they will thus be gen¬ 
erated regardless of the envelope mass (although the enve¬ 
lope can influence their geometric development). Since line 
polarization indicates patchy absorption, the comparatively 
high strength of Ha polarization in SN 2008ax could be re¬ 
lated to the compactness and relatively thin H envelope of 
the progenitor, whereby plumes of ®®Ni from the core might 
retain their jet-like morphology with less overlying mass to 
decelerate and mix them, giving rise to the higher degree of 
axisymmetry. Alternatively, the relatively fast radial veloci¬ 
ties of line polarization in SN 2008ax might also imply that 
a ®®Ni plume or jet was more highly aligned with our line 
of sight and thus preferentially obscured the more forward- 


scattered (less polarized) photons near the apparent centre 
of the photosphere, giving rise to a higher degree of line 
polarization. Of course, such viewing-angle effects could po¬ 
tentially blur the evidence for any intrinsic geometric simi¬ 
larities that SNellb as a class might share, whether they be 
tidal deformation effects, aspherical plumes of radioactive 
®®Ni, or something else (e.g., the potential effect of oblique 
shock breakout discussed in §4.1.2). 

Comparison of different SNellb is further complicated 
by the potential presence of an additional component of con¬ 
tinuum polarization from scattering of optical photons from 
dust in pre-existing CSM. Indeed, SN 1993J exhibited a rel¬ 
atively large UV excess on day 19, attributable to a substan¬ 
tial amount of CSM, while SN 201 Idh exhibited insignificant 
UV excess at a similar phase (Ben-Ami et al. 2015). Scat¬ 
tering by dust might thus be the most plausible explanation 
for the relatively large and long-lasting continuum polariza¬ 
tion in SN 1993J, a possibility that was discussed previously 
by Hoflich et al. (1996) and Tran et al. (1997). The more 
recent SN 2013df also appears to have had a substantial UV 
excess from CSM (Ben-Ami et al. 2015), and it will be inter¬ 
esting to see whether the spectropolarimetric data of that 
event share any noteworthy similarities or differences with 
SN 1993J and SN2011dh that might be attributable to dif¬ 
ferent amounts of CSM in each case. If hypothetically added 
to SN 2011dh, another component of continuum polarization 
from dust scattering could have potentially wiped out the 
approximate symmetry between the position angles of the 
continuum polarization at early times and line polarization 
at later times. Complementary UV data could thus provide 
important constraints for the interpretation of spectropo¬ 
larimetry, helping one to anticipate the potential effects of 
CSM dust. 


4-2.1 Insight from Cas A 

The scenario of an aspherical SN Hb explosion influenced by 
jets or plumes of enhanced radioactive heating is consistent 
with the detailed structure of the Cas A supernova rem¬ 
nant (SNR) and the explosion itself. Spectra of Cas A light 
echoes have unambiguously shown that this was a SNIIb 
(Krause et al. 2008), and the substantial differences in the 
line strength and velocity from multiple echoes that mir¬ 
rored the SN from different viewing angles is direct evidence 
for an aspherical explosion (Rest et al. 2011). Interestingly, 
the SNR also exhibits a highly aspherical distribution of O 
and Fe knots, which are thought to have been produced by 
jets or plumes of radioactive ®®Ni that penetrated the outer 
layers of the SN at early times, or from the growth of RT in¬ 
stabilities (Milisavljevic & Fesen 2013, 2015). In addition to 
our spectropolarimetry, the strong multi-peaked [O l] emis¬ 
sion detected from SN2011dh at late times by Shivvers et 
al. (2013) provides independent evidence that an aspherical 
0-rich SNR like Cas A could be developing. Chornock et al. 
(2011) drew similar parallels between Cas A and SN 2008ax. 


5 SUMMARY AND CONCLUDING REMARKS 

We present seven epochs of spectropolarimetry of the Type 
Hb SN2011dh in M51. The observations cover 86 days of 
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evolution, from the early H-rich stage through the nebu¬ 
lar phase. This represents the most complete spectropolari- 
metric coverage obtained for a SNIIb since SN 1993J. The 
late-phase polarization at the central wavelengths of nebu¬ 
lar features and the q-u loops of Call IR lines at all phases 
indicate a low value of total ISP (< 0.1%), and we have thus 
made no attempt to remove it from the data. At each epoch, 
considerable enhanced polarization is detected across line 
features, while substantial continuum polarization is hrmly 
detected on days 9 and 14. Several physical scenarios could 
potentially explain the spectropolarimetric evolution, sum¬ 
marised as follows. 

The global asphericity of the early-time photosphere in¬ 
ferred from the continuum polarization could be the imprint 
of tidal deformation of the YSG progenitor by the purported 
massive companion; this hypothesis is facilitated by the re¬ 
cent detection of a hot star at the explosion site by Folatelli 
et al. (2014), while substantial geometric deformation of the 
progenitor envelope is consistent with physical parameters 
of the binary suggested by Benvenuto et al. (2013). This hy¬ 
pothesis does not, however, provide a clear explanation for 
the approximate symmetry between the early continuum po¬ 
larization and the subsequent Ha and Hel line polarization 
that we observe (see §4.1.3). 

Alternatively, the continuum and line polarization could 
be the imprint of clumpy excitation of the envelope by fast¬ 
rising plumes of ®®Ni from the core. This interpretation does 
has the advantage of explaining the apparent symmetry be¬ 
tween the early-time continuum polarization and the subse¬ 
quent Ha and Hel line polarization. We currently lean to¬ 
ward this scenario (at least for explaining the line polariza¬ 
tion) , as jets and plumes of ®®Ni are well motivated by hydro- 
dynamic models (Iwamoto et al. 1997; Hammer et al. 2010; 
Wongwathanarat et al. 2010) and have commonly been in¬ 
voked to explain the spectropolarimetric and nebular-phase 
data of other CC SNe, and in the interpretation of detailed 
structure in the Gas A SNR. 

For OI polarization, the independent geometry and 
kinematics suggest that the inferred dumpiness of the ab¬ 
sorption could be the result of RT fluid instabilities that de¬ 
velop near the He/C-fO interface, which is consistent with 
predictions from hydrodynamic models of SNellb (Iwamoto 
et al. 1997). For Ca ll IR polarization, the independent geom¬ 
etry and relatively high-RV edge of the polarized line profile 
at early times potentially suggests an origin in the outermost 
regions of the outflow, perhaps in the clumpy swept-up wind 
of the progenitor. 

Another potential explanation for the early-time con¬ 
tinuum polarization is optical scattering of SN photons by 
dusty CSM. Indeed, the presence of extended CSM is in¬ 
dicated by X-ray observations of SN 2011dh out to day 
~ 500, consistent with ongoing interaction between the re¬ 
verse shock and CSM produced by a steady-state progenitor 
wind (Maeda et al. 2014). However, if dust scattering is re¬ 
sponsible for the continuum polarization, then it is puzzling 
why SN 201 Idh has such a meager UV excess compared with 
SN 1993J (Ben-Ami et al. 2015), since a mass of CSM capa¬ 
ble of scattering the radiation and producing strong contin¬ 
uum polarization should have produced a non-negligible UV 
excess. Furthermore, the fact that our shortest-wavelength 
coverage of the polarization on days 14 and 58 does not ex¬ 
hibit an systematic increase in the magnitude of polarization 


appears to be inconsistent with scattering of SN photons off 
of dusty CSM. Such scattering by dust could have been mit¬ 
igated if grains in the CSM were effectively destroyed by 
the radiation from the SN during the early phases of the 
explosion. 

An aspherical explosion stemming from oblique shock 
breakout is another possible source of continuum polariza¬ 
tion in SN 2011dh. Theoretical calculations indicate that 
this phenomenon might be relatively common for stripped- 
envelope SNe (Matzner et al. 2013); perhaps this could ex¬ 
plain the commonality of polarized flux from such SNe. 

Although the above interpretations of the spectropo- 
larimetry are plausible and guided by theoretical mod¬ 
els, more sophisticated quantitative analysis is required to 
further discriminate between them and explore potential 
combinations of their effects. Specifically, the construction 
of more detailed and robust hydrodynamic and electron¬ 
scattering models that account for the geometries we have 
discussed will enable more accurate interpretations of the 
emergent polarization. 

SNellb are not as observationally homogeneous as pre¬ 
viously thought. However, some interesting and potentially 
illuminating patterns appear to be emerging for this class. 
The frequent detection of substantial polarization demon¬ 
strates that asphericity is common, and thus that the as¬ 
sumption of spherical symmetry is dubious; this is probably 
the reason why SN 201 Idh was classified as a cllb subtype 
on the basis of its early-time luminosity and radio emis¬ 
sion (Arcavi et al. 2011), even though the explosion clearly 
stemmed from an extended YSC (Van Dyk et al. 2013). 
Additional observations of the cllb and ellb subclasses un¬ 
dergoing the transition from H-rich to He-rich phases, and 
complementary UV data indicating the potential for CSM 
interaction and dust scattering, will be necessary to more 
firmly establish any systematic trends in their spectropo¬ 
larimetric properties and to determine the physical reasons 
for their differences. As it stands, explosion energetics, ra¬ 
dioactive inhomogeneities, the mass and potential deformity 
of stellar envelope, optical scattering by CSM dust, the size 
and distribution of absorbing clumps relative to the pro¬ 
jected photosphere, and viewing-angle effects, are all likely 
to be important factors in the increasing spectropolarimetric 
variety of SNellb. 
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